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UNIDIRECTIONAL EVOLUTION EQUATIONS OF DIFFUSION TYPE 


GORO AKAGI AND MASATO KIMURA 


Abstract. This paper is concerned with the uniqueness, existence, comparison principle and 
long-time behavior of solutions to the initial-boundary value problem for a unidirectional dif¬ 
fusion equation. The unidirectional evolution often appears in Damage Mechanics due to the 
strong irreversibility of crack propagation or damage evolution. The existence of solutions is 
proved in an L^-framework by introducing a peculiar discretization of the unidirectional diffu¬ 
sion equation by means of variational inequities of obstacle type and by developing a regularity 
theory for such variational inequalities. The novel discretization argument will be also applied 
to prove the comparison principle as well as to investigate the long-time behavior of solutions. 


1. Introduction 

In Damage Mechanics, the unidirectional evolution is a significant feature of crack propagation 
models. Indeed, crack propagation is an irreversible phenomena, and particularly, cracks in 
a specimen or the damage of a material (e.g., microcracks which break or weaken bonds of 
microstructures) cannot disappear nor decrease. Hence if one introduces a phase parameter 
(see [Ml Ea SSI EH Eg ) or an internal variable (see mm) which describes the crack growth 
or the damage accumulation, they are usually supposed to be unidirectional, i.e., nondecreasing 
or nonincreasing. Such unidirectional evolution processes are often described by PDEs involving 
the positive-part function, s i-A (s)+ := s V 0 = max{s, 0} for s G R. Moreover, in a phase field 
approach, the evolution of a phase parameter is governed by the gradient flow of an appropriate 
energy functional subject to some unidirectional constraint. 

In order to find out mathematical features of such unidirectional evolutions and to develop 
mathematical devices to analyze them, in this paper, we shall treat, as a simplest case, the 
evolution of u = u{x, t) governed by the (fully nonlinear) PDE, 

9tM=(Au + /)^, for X G fl, t > 0, (1) 

where 17 is a bounded Lipschitz domain of R" with n G N, dtu = du/dt, A stands for the n- 

dimensional Laplacian, / = f{x,t) is a given function and (s)+ := s VO for s G M. More precisely, 

the main purpose of this paper is to prove the uniqueness, existence and comparison principle 
of strong solutions u = u(x, t) of the initial-boundary value problem for ([T]) and to reveal the 
asymptotic behavior of u = u{x, t) as t —>■ oo. 

Solutions of o entail unidirectional nature, more precisely, the non-decrease of m = u{x,t) in 
t, since the right-hand side of m is non-negative due to the presence of the positive part function. 
There appear various kinds of unidirectional evolutions in natural sciences and engineering fields 
(see, e.g., [21]). In particular, a phase field model for crack propagation in an elastic material was 
proposed with the aid of a unidirectional gradient flow by the second author and his collaborator 
in [MIIIZIIIS] (see also [BJED])- Equation ([T]) can be regarded as a simplified equation of their 
model. 

In mathematical points of view, © is classified as a fully nonlinear PDE, which is not fit for 
energy methods in general; however, by taking a (multi-valued) inverse function of the positive 
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part function (■ )+, o can be formulated as a sort of doubly nonlinear evolution equations, 

dtu + a{dtu) — Au ^ f a.e. in x (0, oo), (2) 

where a is a (multi-valued) maximal monotone function in R given by a(0) = (—oo, 0] and 
q;(s) = {0} for any s > 0 with the domain D{a) = [0, oo) (see Remark [2?3l and Section [4] below for 
more details). Equation JS]) is fitter for energy methods and monotone techniques. On the other 
hand, in view of the L^-theory of evolution equations, two operators v i—>• a{v{-)) and v !—>■ —Av 
(defined for v € are unbounded in and hence, it is more delicate to establish a 

priori estimates for proving the existence of strong solutions, as compared with standard equations 
without unidirectional constraints, e.g., the classical and nonlinear diffusion equations. 

The nonlinear PDE @ may fall within the frame of abstract doubly nonlinear evolution 
equations in a Hilbert space H of the form 

d'^{dtu{t)) -I- d^{u{t)) 9 f{t) in i7, 0 < t < T, (3) 

where ST and 94) denote the subdifferential operator of functionals H ^ (—oo, oo] and 
d* : H > (— 00 , 00 ], respectively. In a thermodynamic approach to continuum mechanics, \E' and 
d> are often referred to as a dissipation functional and an energy functional, respectively. To 
reduce (I2|) into the form Q, we set u{t) := u{-,t) and particularly choose 

H = L^{n), 4'(r;) = i f |ul^da:-b/[ .>o](?^), ^’(?^) = ^ f jVupda; for v G H, 

where /[ . > 0 ] is the indicator function over the set {?; G L'^(yi ): u > 0 a.e. in H}. Then we note 
that both subdifferentials 94' and 94> are unbounded in H. 

Let us briefly review the previous studies on abstract doubly nonlinear evolution equations 
such as Barbu m proved the existence of solutions for Q with two unbounded operators 
94' and 94> by using the elliptic-in-time regularization and by imposing the differentiability (in 
t) of /. This result was generalized by Arai [7], Senba [ID] and so on. In these papers, the term 
94'(9tit(t)) is estimated by differentiating the equation and by testing it with 94'(9tit(t)). There¬ 
fore the differentiability of / (more precisely, / G (0, T; iJ) by |7|) is essentially required, and 
moreover, some strong monotonicity condition (i.e., the so-called 94'-monotonicity) is also im¬ 
posed on 94>. Similar methods of establishing a priori estimates are also used in individual studies 
on irreversible phase transition models (see, e.g., [14]). On the other hand. Colli and Visintin es¬ 
tablished an alternative approach to ([3]) in [18] , where 94' is supposed to be bounded and coercive 
with linear growth instead of assuming the regularity assumption on / and the 94'-monotonicity 
of 94’ (see also m)- Their framework would be more flexible in view of applications to nonlinear 
PDEs and has been extensively applied to various types of doubly nonlinear problems. Moreover, 
their framework has been generalized in many directions, e.g., perturbation problems, long-time 
behaviors (see, e.g., [UIMIISS], [DU Sect. 11], [HKMIIMII^[111 1113). However, due to the 
unboundedness of 94', ([D]) seems to be beyond the scope of the latter approach. On the other 
hand, the former approach due to Barbu and Arai is applicable to ([D]); however, the regularity 
assumption / G W^’^{0,T; H) has to be assumed for proving the existence of solutions. Aso et 
al. [HE] also treated an irreversible phase transition system in a different fashion; however, the 
regularity condition / G H) is also assumed there. 

In this paper, we present a novel approach to © ( or equivalently ©) by introducing a peculiar 
discretization for 0 by means of variational inequalities of obstacle type. Moreover, by developing 
a regularity theory of such variational inequalities, we shall establish new a priori estimates for 
(or ©) without assuming the differentiability (in t) of /. Such a relaxation of the regularity 
assumption on / may bring some advantage to consider perturbation problems, which will be 
left for a forthcoming paper. Moreover, the novel discretization method will be also applied 
to investigate the long-time behavior of solutions as well as to prove a comparison theorem. 
In particular, we shall provide uniform (in t) estimates for solutions by employing the peculiar 
discretization scheme. Eurthermore, some variational inequality of obstacle type will play a 
crucial role in asymptotic analysis; indeed, it will turn out that every solution will converge to 
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the unique solution z = z{x) of a variational inequality of obstacle type involving the initial data 
as an obstacle function from below under suitable assumptions. Here it is worth mentioning that 
the limit z of the solution u = u{x,t) depends on its initial data uq; indeed, one can construct 
different limits of solutions for different initial data. On the other hand, the w-limit set of each 
solution must be singleton. 

The organization of this paper is as follows. The initial-boundary value problem for m and 
basic notions (e.g., strong solution) and assumptions are formulated, and main results of this 
paper will be stated in the next section. Section [3] is devoted to establishing a regularity theory 
as well as to verifying a comparison theorem for variational inequalities of obstacle type. In 
Section 01 we discuss a reduction of m to an evolution equation of doubly nonlinear type in 
L^(r2) and prove the uniqueness of solutions for the initial-boundary value problem. In Section 0 
we carry out the backward-Euler time-discretization of © and construct a strong solution of 
(HI) by proposing a new a priori estimate based on the regularity theory developed in Section 01 
A comparison theorem for is also proved here. The long-time behavior of solutions will be 
investigated in Section 01 In the last section, we shall discuss other equivalent formulations of 
solutions for ©• 

Notation. For each normed space N, we denote by N' the dual space of N with the duality 
pairing {g, v) n := n' {g, v)n = g{v) for u G N and g G N'. For Banach spaces U and W, the set of 
all bounded linear operators from U to W is denoted by B{U, W). Moreover, the set of all linear 
topological isomorphisms from U to W is denoted by Isom([/, W), that is, A G Isom(C7, V) means 
that A is bijective from U to W, A G B{U, W) and A~^ G B(W, U). Furthermore, Ti.^ stands for 
the /c-dimensional Hausdorff measure in R" for fc = 1, 2,.. ., n. We also write aV b = max{a, b} 
and aAb = min{a, b} for a, 6 G R. Finally, C denotes a non-negative constant independent of the 
elements of the corresponding space and set and may vary from line to line. 

2. Main results 

Let H be a bounded Lipschitz domain in R" with n G N. Let T be the boundary of such 
that r is the disjoint union of two subsets Td and Tn, that is, 

rDurN = r, rDnrN = 0. 

Moreover, assure that T^i is (relatively) open in L. Let v denote the outward-pointing unit 
normal vector on T. One of these two subsets may be empty. In such a case, the other set 
coincides with the whole of L. Main results of the present paper are concerned with the following 
initial-boundary value problem for a unidirectional evolution equation of diffusion type. 


dtu = (Am + /)_^ 

in Q := fi X (0, oo), 

(4) 

u = 0 

on Td X (0, oo), 

(5) 

duU = 0 

on Tn X (0, oo). 

(6) 

II 

o 

II 

e 

o 

in B, 

(7) 


where dt = d/dt, f = f{x,t) and ug = uo(x) are given functions of class Lf^^{[0,oo)-, and 

L^(H), respectively, and d^u := 7 o(Vit) • u denotes the normal derivative of u. Moreover, (•)+ 
stands for the positive part function, i.e., (s)+ := s V 0 for s G R. If Ld (resp., Tn) is empty, the 
corresponding boundary condition (0]) (resp., (01)) is ignored. 

Remark 2.1. By change of variable, one can reduce another unidirectional diffusion equation, 

dtU={Au + f)_ in Q, (8) 

where (s)_ := s A 0 for s G R, to (01). Indeed, set v := —u and g := —/. Then ([8]) is transformed 
to 

-dtv = {- Av-g)_ = -{Av + g)^, 
whence v solves (0]) with / replaced by g. 
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Let US start with defining strong solutions of O-©. To this end, we set up notation. Let 
7 o G iL^/^(r)) denote the trace operator defined on iL^(O) (throughout the paper, we 

may omit 70 if no confusion can arise). Moreover, define 

V := {v € 70 U = 0 'H"“^-a.e. on Td}, 

X := {v G H^{n): 7 o(Vu) ■ z/= 0 ■H"-i-a.e. on Tn} 

equipped with the induced norms and inner products, i.e., IHly = H'lli/qn) and(-,-)v = 
for V] II ■ ||a = II • ||ff 2 (n) and (•, ■)x = (•, ■)H'^{n) for X. Then V and X are closed subspaces of 
and iL^(n), respectively; hence, they are Hilbert spaces. If either Td or Tn is empty, the 
corresponding boundary condition specified in the definition of H or X above is ignored. 

We are concerned with strong solutions of ©-© defined by 

Definition 2.2 (Strong solution). For T > 0, a function u G C([0, T]; L^(H)) is called a strong 
solution of 0 -© on [ 0 ,T], if the following three conditions are satisfied: 

(i) u G Wi’ 2 ( 0 , T; n L^{0, T;Xn V), 

(ii) the equation dtu = (Au + /)+ holds a.e. in Qt := H x (0,T), 

(iii) the initial eondition u|t=o = uq holds a.e. in H. 

A function u G CdO, 00 ); L^(H)) is called a strong solution of ©-© on [0, 00 ), if for any T > 0, 
the restriction of u onto [0,T] is a strong solution of ©-© on [0,T]. 

Remark 2.3. One can further derive that u G C'([0,r];H) from (i) by employing a chain-rule 
for convex functionals. See Lemma 2^ in (2] below for more details. 


We are now in position to state main results, whose proofs will be given in later sections. We 
begin with the uniqueness of solutions. 

Theorem 2.4 (Uniqueness). Let T > 0, uq € V and f G L^{Qt). Then the strong solution of 
0-0 on [0, T] is unique. 


To state our existence result, we shall introduce some assumptions for the domain and the 
boundary Td, Tn- For A G K, we define a mapping A\ G B{V, V') by 


{A\u,v)v = / (Vit • Vu -I- Auu) dx for u,vGV. 
Jn 


(9) 


It is well known that A\ G Isom(y, V) holds if A > 0. Hence one can define u = A^ g for 
g G as the unique solution u of the elliptic problem in a weak form, 

/ (Vit • Vx -I- \uv) dx = / gv dx for all v G V, 

Jn Jn 

(i.e., A\u = g in V'). Then we assume that 

Af^g G for all g G (10) 


Condition (HOI) is often called an elliptie regularity condition and deeply related to the geometry 
of the domain and boundary conditions. Indeed, it holds true for smooth domains with a single 
boundary condition (i.e., rN = 0orrD = 0). However, it is more delicate to consider the 
validity of HOD for situations with nonsmooth domains or mixed boundary conditions. On the 
other hand, in order to take account of physical backgrounds of crack growth models and their 
numerical simulations, the regularity of the boundary may be at most Lipschitz continuous, and 
mixed boundary conditions seem to be natural as well. We shall give conditions equivalent to 
in Proposition 13.61 below. 


Remark 2.5. Let us exhibit a couple of examples of H, Tp and Tn for which the condition (flUl) 
is satisfied. 


(i) If Tp n Tn = 0 and L is of class then (ITUl) is satisfied (see Theorem 2.2.2.3 and 

Theorem 2.2.2.5 of [25)1. 


UNIDIRECTIONAL EVOLUTION EQUATIONS OF DIFFUSION TYPE 


5 


(ii) Let il be convex. If Ld = L or Ld = 0, then (flOl) is satisfied (see Theorem 3.2.1.2 and 
Theorem 3.2.1.3 of [25)'). 

(iii) If n = 1 or if 17 is a rectangle in and Ld is a union of some of four edges of T, then 

is satisfied. Indeed, since the weak solution u = g can be extended to an open 
neighborhood of 17 by reflection, it follows that u = A^^g € 

Our existence result reads. 

Theorem 2.6 (Existence). We suppose that the condition (1101) holds true. LetT > 0, uq G XCiV 
and f G L^{Qt) be given. If there exists f* G satisfying 

fix,t)<f*{x) a.e.inQr, (11) 

then there exists a strong solution u = u(x,t) to the problem ©-([T]) on [0,T]. 

Remark 2.7 (Assumptions on /). Condition (fTT|l is weaker than / G L°°{Qt) or f G T ; T^(17)) 

(cf. [7]). Indeed, if / € 1T^’^(0, T; L^(17)), then f*{x) := f{x, 0)+ \dtf(x, t)| dt belongs to L^(17) 
and satisfies (ITTI) . On the other hand, (fTTI) is stronger than /+ := / V 0 € L°°(0, T; L^(17)). In 
fact, (fTT|) yields /+ G T°“(0, T; L^(17)). However, even if /+ G L°“(0, T; L^(17)), (|TT|l might not 
hold true. One may easily find a counterexample, e.g., f{x,t) = |a; — t|““, 17 = (0,1), T = 1 and 
0 < a < 1/2. 

Theorem 12.61 will be proved in Section [S] Our method of proof relies on the backward-Euler 
time-discretization of (|T|) . The discretized equation will be rewritten as some variational inequal¬ 
ities of obstacle type. Moreover, developing a regularity theory for such variational inequalities, 
we shall obtain a new a priori estimate for discretized solutions. As a result, we shall construct 
a strong solution without assuming the differentiability of / in t, which is usually assumed in 
standard frameworks of doubly nonlinear evolution equations such as dll]) (see, e.g., [niiiiiiQ]). 

Moreover, such a basic strategy of proving the existence result will be also applied to prove 
the following comparison theorem for strong solutions of 0-0: 

Theorem 2.8 (Comparison principle). Let T > 0 and suppose that (11011 is satisfied. For each 
i = 1,2, let Uq G X nV and /* G L^{Qt) be such that there exists f* G L^{Q) satisfying 

/*(x,7) < f*ix) a.e. in Qt. 

For i = 1,2, let rt® = u'^(x,t) be the unique strong solution of 0-0 with ug = Uq and f = P 
on [0,T]. If Uq < Uq a.e. in 17 and p < p a.e. in Qt, then P < a.e. in Qt. 

The peculiar discretization argument will be also employed to investigate the long-time behav¬ 
ior of solutions for 0-0. Furthermore, the comparison theorem stated above will play a crucial 
role to identify the limit of each solution u = u(x, t) as t ^ oo. 

Theorem 2.9 (Convergence of solutions as t ^ oo). Let ug G X r\V and assume that (flOll holds 
and that 

(HI) n^-pTo) >0; 

(H2) there exists a function /oo G LfiLl) such that f — foa belongs to T^(0, oo; L^(17)); 

(H3) / G L°°(0, oo; T^(17)), and (ITT]) is satisfied. 

Then the unique solution u = u{x, t) of 0-0 on [0, oo) converges to a function z = z{x) G XDV 
strongly in V as t ^ oo. Moreover, the limit z satisfies 

z > Uq and — Az > /oo a.e. in 17. 

In addition, if f{x,t) < fooix) for a.e. {x,t) G Q, then the limit z coincides with the unique 
solution z G X nV of the following variational inequality. 

z G Kq{uq) -.= {v GV : V >Uq a.e. in 17}, 


(VI)(uo,/oo) 


Vz-X{v — z)dx> / foc{u — z)dx for all v G Kq{uq). 
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Remark 2.10. Assumption (HI) is essentially required to ensure the convergence of the solution 
u = u{x,t) as t ^ oo. Indeed, suppose that Fd = 0 (i.e., Fn = F) and set uo{x) = 1 and 
f{x,t) = 1. The unique strong solution of (H])-© is given by 

u{x,t) = l + t for {x,t) € Q, 

and then, u{x, t) is divergent to oo at each a; S H as t —>■ oo. 

3. Regularity theory for variational inequalities of obstacle type 

In this section, based on the approach of Gustafsson [2^, we revisit a regularity theory for 
variational inequalities of obstacle type. In classical literature on variational inequalities of obsta¬ 
cle type, the IT^’P(H)-regularity of solutions is often obtained by using a penalization technique 
(see, e.g., [551 Chap. 4], [55]). On the other hand, Gustafsson [55] gave a simpler alternative proof 
by introducing an auxiliary variational inequality and by proving the coincidence of solutions for 
both problems. Let us also remark that, in previous results, it is assumed that LF^’^(f2) C (7(0) 
(namely, p > njl') in order to utilize the classical maximum principle for linear elliptic equations. 

We shall establish a VF^’^(0)-regularity result for variational inequalities of obstacle type 
equipped with a mixed boundary condition by properly modifying the argument of Gustafs¬ 
son [5^. It is noteworthy that we do not assume that VF^’P(O) C (7(0), i.e., p > n/2, as we 
employ Stampacchia’s truncation technique instead of the classical maximum principle. 

First, we shall set up notation. For tr > 0, let A := A^, £ i?(I/, V) be defined as in I©. Define 
a symmetric bilinear form a(-, •) : H x H R associated with A by 

a{u,v) := {Au,v)v = / {S/u-Vv + auv) dx for u,v£V. (12) 

Jn 

Throughout this section, we assume that 

(T>0 if H”“^(Fd) = 0. (13) 

Under the condition (1151) . (by the Poincare inequality for the case that cr = 0 and H"“^(Fd) > 0), 
a(-,-) turns out to be coercive on U x U. Hence A is invertible, and A~^ belongs to B{V',V). 
Let / £ U' and ip €V and define closed convex subsets Kq, Ki of U by 

Kq := {v G V: V > Ip a.e. in O}, (14) 

Ki := {v G V: Av > f in V'}. 

Here the inequality Av > f in V' means that {Av — f,<p)v > 0 for ali ip G V satisfying p > 0 
a.e. in D. We also define functionals J,J on V by 

J{v) :=^a{v,v) - {f,v)v for i; £ U, (15) 

J{v) :=^a{v,v) - {f,v)v ioi v G V, 

where / := Aip G V. Then the following equivalence holds true: 

Proposition 3.1. Suppose p3l) and let f G V' and ip G V. Then the following five conditions 
for u G V are equivalent to each other: 

(a) u G Kq, J{u) < J(v) for all v G Kq, 

(b) u G Kq, a{u, V — u) > {f,v — u)v for all v G Kq, 

(c) u G Kq n Ki, {Au — f,u — ip)v = 0, 

(d) u £ Ki, a{u, V — u) > {f,v — u)v for all v G Ki, 

(e) u G Ki, J{u) < J{v) for all v G Ki. 

Moreover, there exists a unique element u G V satisfying all the conditions. 
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Proof. Since J is a coercive, continuous, strictly convex fnnctional on the closed convex set Kq, 
J admits a unique minimizer u over Kq. Hence u satisfies (a). 

We shall prove the equivalence of the conditions (a)-(e). It is well known (see, e.g., [5S]) that 

(a) (b) and (d) (e). So, let us here start with showing that (b) => (c). The condition (b) is 

equivalently rewritten by 

u G Kq, {Au — f,v — u)v > 0 for all v G Kq. (16) 

For any (f G V with (p > 0 a.e. in H, substituting v = u+p G Kq to (fTHll . we have {Au — f, (p)v > 0, 
which yields that u G Ki. On the other hand, substitute v = ip G Kq and v = 2u — tp G Kq to 
Then one can obtain {Au — f,u — 'ip)v = 0. Hence (c) holds. 

To prove the inverse relation, (c) ^ (b), let u satisfy (c). For any v G Kq, we see that 

a{u, V - u) - {f,v - u)v = {Au - f,v - ip)v - {Au - f,u- ip)v. 

Since the first term of the right-hand side is nonnegative (by u G Ki and v G Kq) and the second 
term vanishes (by the equation of (c)), the condition (b) follows. 

One may prove the equivalence between (c) and (d) in a similar fashion to the above. Firstly, 
suppose that u satisfies (c). For any v G Ki, one finds that 

a{u, V — u) — {f,v — u)v = {Av — Au, u)v — {Av — Au, 'ip)v 

= {Av - f,u- ip)v - {Au - f,u- ip)v 

Here we used the fact that {Aw,z)v = {Az,w)v for all w,z G V. Noting that the right-hand 
side is non-negative by (c) and the fact that u G Kq and v G Ki, one can get (d). To check the 
inverse relation, we also rewrite (d) as 

u G Ki, {Av — Au,u — tp)v > 0 for all v G Ki. (17) 

For any p G LF'{VI) with p > Q a.e. in 0, substituting v = u + A~"^p G Ki to (flTll . we have 
{p, u — Pi 0, which along with the arbitrariness of (/? > 0 implies u G Kq. Moreover, let us 

also substitute v = A~^f G Ki and v = 2u — A~^f G Ki in (flTl) . Then we obtain {Au — f,u — 
ip)v = 0, whence (c) follows. Consequently, all the conditions (a)-(e) are equivalent. □ 


In the rest of this section, let p G M satisfy 

1 < p < oo, p > 


2n 


(18) 


Since the Holder conjugate q := p/(p — 1) of p satisfies q < 2* := 2nf{n — 2) if n > 3, and H 
is a Lipschitz domain, by Sobolev’s embedding theorem, the continuous embeddings V ^ L‘^{fl) 
and LP{n) = (T'?(H))' ^ V hold true. We also note that W‘^'P{Vt) is continuously embedded in 
H^{Vt) by (IT^ . 

We further suppose that 

fGLP{n), iPgV, AxPGLP{n) (19) 

and introduce a closed (in V) convex set K 2 given by 

K 2 ■= {v G V : f < Av < / V / in V'} C Ki. 


A main result of this section is stated as follows: 


Theorem 3.2. Suppose that (11311 . (1181) and (1191) are satisfied. Then each of the following condi¬ 
tions (f)-(h) is equivalent to the conditions (a)-(e) of Provosition 

(f) u G K 2 , J{u) < J{v) for all v G K 2 , 

(g) u G K 2 , a{u, V — u) > {f,v — u)v for all v G K 2 , 

(h) u G Kq in K 2 , {Au — f){u — Ip) = 0 a.e. in H. 


We first prepare a couple of lemmas, which will be also used in later sections. 
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Lemma 3.3. Let w G and set w+(x) := (wix))+ for x € id. Then w+ G H^(il) and 

7o(w+) = (7oR')+ 'hr~^-a.e. on V. 


Proof. Set := {cc G w{x) > 0} and recall Theorem A.l.of | 
and 


to observe that w+ G H^{il) 


w+ = 


w a.e. in f2+, 
0 a.e. in \ 


Vw+ = 


Vw a.e. in f2+, 

0 a.e. in id \ n+. 


( 20 ) 


Set W := (7(0) n H^{Q). Since W is dense in iJ^(O), there exists a sequence {wn} in W such 
that Wn ^ w strongly in i7^(0) as n oo. Noting that 

||(ui„)+- ^+11^2(0) < - w||i2(Q) 0 as n-)> oo, ( 21 ) 

we observe that {wn)+ w+ strongly in L^(0). Applying (1^ to w;„, we also have 

II ('R’n)+= ll('fi^n)+llL2(f2) + || V(r(;„)+||^2 (q) < ll'^n|li2(Q) + II VlC„||^2(Q) = ||R'n||//i(n)' 

Since {iCn} is bounded in i7^(0), so is {(w„)+}. Hence, one can extract a (non-relabeled) 
subsequence of {n} such that {wn)+ —>• w+ weakly in i7^(0) Again from (l20l) . we have 

||(t«n)+||//i(n) = ll(^"n)+||i2(Q) + (V(ii;„)+, Vw„)L2(n) 

^ Il^+lli2(n) + (Vw+, Vu;)L2(f2) = ||R'+||^i(n): 

which together with the uniform convexity of iJ^(O) also implies that {wn)+ —^ '!c+ strongly in 
H\n). Since 70 G H(iJi(0),iJi/2(r)), we particularly deduce that yoWn ^ 'jow and 7o(rCn)+ —?• 
7 ow+ strongly in A^(r). As in (l2T]l . one can verify that (7oWn)+ —t (7o'R')+ strongly in A^(r). 
On the other hand, since Wn G (7(0), it is clear that 

7o(w„)+ = {loWn)+ H”“^-a.e. on T. (22) 


Passing to the limit as n —>■ oo in (|22ll . we conclude that 7 orc+ = ( 7 orc)+ ^-a.e. on P. 

Lemma 3.4. If vi,V 2 G V, then viV V 2 & V and vi A V 2 € V. 


□ 


Proof. Applying Lemma [3.31 to w = ±(ui — V 2 ), we find that {vi — V 2 )+ and {v 2 — wi)+ belong 
to V. Hence we obtain z;i V U 2 G H and z;i A U 2 G V, since 'Ui V n 2 = U 2 + (wi — ^ 2 )+ and 
til A r;2 = ^^2 — {v2 — vi)+. □ 


Let us move on to a proof of Theorem 13.21 

Proof of Theorem, \S.^ It is obvious that (f) <(=> (g). As in Proposition 13.11 one can uniquely 
choose u € K 2 which satisfies (f) and (g). Next, we shall prove the equivalence between (a)-(e) 
and (f), (g). Let ui be the unique element of V satisfying (a)-(e) and let U 2 be the unique element 
of V satisfying (f) and (g). 

We claim that Ui = U 2 . Indeed, note that Au 2 G L^{Tl) by U 2 € K^. Set w := U 2 —ip €V and 
h := Aw G LP{il). Since U 2 satisfies (g), it follows that 


0 < a{u 2 ,v — U 2 ) -{f,v — U 2 )v = {Au 2 — Alp, v - U 2 )v 

= {Av — Au 2 , w)v for all v G K 2 . 

We set a measurable set 

N -.= {x & H: w{x) < 0}. 

Define g G L’^iLl) by 


( 23 ) 


( ^ / /(^) V f{x) if X € N, 

ifxen\jv. 
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Then by definition one can observe that A G if 2 - Hence substituting v = A to (IMl) . we 
have 

0 < {g - Au 2 ,w)v = j {{f y f) - Au-^wAx. 

Since {f V f) — Au 2 > 0 and w < 0 a.e. in N, one can derive the relation Au 2 = f ^ f a.e. in N. 
It follows that 

h = Au 2 - = (/ V /)-/ = (/-/)+> 0 a.e. in N. 


On the other hand, we recall that w G V solves the equation Aw = h in H', that is, 


a{w,v) = (h,v)v for all v G V. 

We define W- := w AO. Then by Lemma 15^ and Theorem A.l. of [35], we have 

I w in TV, I Vw a.e. in N, 

W- G V, W- = < , vw- = < 

|0 in 0 \ TV, |0 a.e. in 0 \ iV. 

Substituting v = W- into (IMl) . we have 

0 < a{w-,w-) = a{w,w-)^^ {h,w-)v = / hwdx. 

JN 


(24) 


(25) 


Since h > 0 and tc < 0 a.e. in N, the right-hand side of (1351) is non-positive. Hence we deduce 
that a{w-,w-) = 0, which along with the coercivity of a(-, •) implies that W- = 0 (i.e., tc > 0) 
a.e. in fi. Therefore U 2 belongs to ifo- 

Substitute v = A~^f G if 2 to the condition (1331) . Then we obtain 


{Au2 - /,U2 - '0)v < 0. 


(26) 


Moreover, noting that Au2 — / > 0 (by U2 G K2) and U 2 — '0 > 0 (by U2 G ifo), we derive 
{Au 2 — f,U2 — 'ip)v = 0 by (I26p . Hence, u = U2 satisfies the condition (c). By uniqueness, we 
obtain ui = M 2 . Thus we have proved that all the conditions (a)-(g) are equivalent. 

Finally, we note that (h) immediately implies (c), since if 2 C ifi. Conversely, let u satisfy (c). 
Then u belongs to K 2 by (f), and hence, Au G LP(H) and 


0 = {Au — f,u — 'il))v = I {Au — f){u — ip) Ax. 
in 

Thus we obtain {Au — f){u — ip) =0 a.e. in H, since Au > / and u > ip a.e. in H. Therefore (h) 
holds. □ 


Thanks to Theorem 13.21 for each solution u of the variational inequality of obstacle type, 

u S ifo, a{u,v — u) > {f,v — u)v for all m G ifo, (27) 

we have obtained an additional information, u G if 2 . In order to more explicitly clarify the 
feature of the additional information, let us introduce the following assumption: 


A^ ^g G VF^’P(H) for all g G LP{n) in case tr > 0 or p <2*, 

Ap^g GW‘^’f^{^) for all g G LP{^) and pG[2*,p] incase tr = 0 and p > 2*, 


(28) 


where 2* = 2n/(n — 2)+ (here we note that (flOl) is a special case of 
find that 


with p = 2). Then we 


Proposition 3.5. Assume that 


holds. Then 

K 2 C W2’P(H). 
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Proof. Let v € K 2 . In case u > 0 or p < 2* (hence V ^ since Av G LP{il) and v G V, 

one can check that v G U'{Vl). Hence A^v = Av + (1 — (t)v G which along with (E5)) gives 

V G W‘^’P{il). In case ct = 0 and p > 2*, noting that Aiv G as above, we find by (E51) that 

V G IH^’^ (H). By iteration argument along with Sobolev’s embeddings, one can finally conclude 

that V G Therefore we deduce that v G as in the former case. □ 

Condition (j28p can be regarded as elliptic regularity of weak solutions for the elliptic boundary 
value problem, 

—Au + u = f in H, u = 0 on Td, = 0 on Fn, 

and it holds true in many cases, e.g., smooth domains with Fn = 0 or Fd = 0 (see, e.g., [24]). 
However, the validity of (1281) is more delicate, if H is not smooth or mixed boundary conditions 
are imposed. So we here explicitly made the assumption. 

To take account of boundary conditions, we further define a subspace of by 

XP := {v G W^’P{n) : 7 o(Vu) • = 0 H""^-a.e. on Fn}. 

Prior to stating a regularity result for (071) . we prepare a proposition, which provides equivalent 
forms of the assumption (051) . 

Proposition 3.6. Under the assumption (1181) . the operator A\\xpr\V restricted onto XP C\V is 
injective and bounded linear from XPCiV into LP{UL), and it coincides with the operator —A + A, 
where A means the Laplace operator from D{A) = XP D V into LP(VL), that is, the Laplacian 
equipped with the Dirichlet and Neumann boundary conditions onT^ and Fn, respectively, in a 
strong form. 

Moreover, the following conditions for LI and Fd, Fn are equivalent to each other: 

(i) there exists A > 0 such that Aj^^g G W^'P{LI) for all g G LP{LI); 

(ii) for any X > 0, it holds that Af^g G W'^’P{LI) for all g G LP(Lt)', 

(hi) there exists A > 0 such that Af^g G XP for all g G LP{LI), and (—A + X) G Isom(A^’ D 
V,LPin)); 

(iv) for any X > 0, it holds that A^^g G Xp for all g G LP{LI), and (—A + A) € Isom(AP n 
v,LP{n)). 

Proof. Denote Bx := Ax\xpnv for A > 0. Then, for u G Xp D V, we observe by Green’s formula, 
which is valid for Lipschitz domains, that 

{Bxu,v)v = / (Vu • Vu + Auu) da: = / (—Ait + Au) u da: for all v GV 
Jo, Jn 

which implies that Bxu = —Ait + Xu and Bx G B{XP D V, LP{LI)), since V (A is dense in 

L'J(r2). Thus we obtain Bx = (—A + A). Moreover, Bx is injective, since so is Ax. 

As for the equivalence of (i)-(iv), we shall show (ii) => (i) => (iii) ^ (iv) => (ii). It is clear 

that (ii) => (i) (and also (iv) => (hi)). We show (i) => (iii). Assume (i), let g G LP{Lt) and set 

u := Af^g G W'^’P{Ll) D V. For all u G H D W^’^{Lt), we have 

/ gv Ax = {Axu,v)v = / (Vu • Vu + Aitu) dx = / (cI„it)i'd'H"“^ + / (—Au + Ait) udx, 

«/o </o 

which implies that u G XP and — Ait(x) + Ait(x) = g(x) for a.e. x G LI. Hence u = Bfj^g. 
Therefore, Bx is surjective from XP D V into LP(LL). By the open mapping theorem, we obtain 
(-A + X) = BxG Isom(AP n V, LP(LI)). 

We next show (iii) (iv). Under the condition (iii), it holds that Bx G Isom(AP (A V, LP{LI)), 
in particular, Bx : XP n U —)> LP{LI) is a Fredholm operator of index zero. For arbitrary ^ > 0, 
we find that = Bx + (/i — A) is a Fredholm operator of index zero from Xp DV to LP{LI) as 
well, since XPUV is compactly embedded in LP{LI). Since B^ is injective (i.e., dimker(B^) = 0), 
we infer that B^ is surjective, and hence also belongs to Isom(A^’ (A V, LP(LI)). Furthermore, 
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for any g € LP{Q.) and /i > 0, the element u = (—A + g) ^g belongs to n V. Hence 

Bfj,u = g, i.e., = g, which implies A~^g = u G XP. Thus (iv) follows. 

It is obvious that (iv) implies (ii) by the definition of XP. Thus we have shown that all the 
conditions (i)-(iv) are equivalent to each other. □ 


Remark 3.7. Under (051) . the assumptions for ip in (1151) is equivalent to ip G XP n V. Indeed, 
lei Ip G V satisfy Aip G LP{Gt). Then as in the proof of Proposition 13.51 one can check that 
Ip G W‘^'P{il). Moreover, by Green’s formula, we hnd that 


/ Aipv dx = {Aip,v)v = / Vip-Vvdx + cr / ipvdx 

J Q, J Vt «/ O 

= f {—Aip + aip)vdx+ ( {duip)v d'H^~^ 

J Q, 

for all u G U n Thus by the arbitrariness of v, we obtain d,^ip = 0 7^"“^-a.e. on Tn, 

whence follows ip G XP. 


Now, we are in position to state a regularity result for (071) as a corollary of Theorem l3.2l This 
corollary will be used for proving Theorem 12.61 in Section [5l 

Corollary 3.8 (Regularity of solutions for variational inequalities of obstacle type). Assume 
that (051) . (051) . / G LP{Q), Ip G XP n V and (1281) are satisfied. Let u gV be the unique element 
satisfying (a)-(h). Then it holds that 

u G XP n Kq, f < Au < f y f a.e. in LI. (29) 


Proof. By Theorem 13.21 the unique element u G V satisfying (a)-(h) belongs to K 2 . Hence, by 
Proposition 13.51 one has u G W‘^’P{LI). Since one can observe that u = Af"^{Au + (1 — cr)u) and 
Au + (1 — a)u G LP{fl), by (1551) along with Proposition 13.61 it holds that u G XP. □ 


We next give a comparison theorem for variational inequalities of obstacle type. 

Theorem 3.9 (Comparison principle for variational inequalities of obstacle type). We suppose 
that da and da are satisfied. For i = 1,2, let fi G LP(LI) and ipi G V be such that Aipi G LP[LI) 

and set Kq := {v G V: v > ipi a.e. in fl}. Let Ui G V be the unique solution of the variational 

inequality: 

Ui G Kq, a{ui,v - Ui) > {fi,v - Ui)v for all v G Kq (30) 

for i = 1,2. If fi < /2 and ipi < ip 2 a.e. in LI, then ui < U 2 a.e. in LI. 

This theorem will be used to prove Theorem 12.81 a comparison theorem for the evolutionary 
problem (jT|), with the aid of the discretization argument. 

To prove Theorem 13.91 we prepare the following lemma. 

Lemma 3.10. We suppose that (1131) . (1181) and (1191) are satisfied. Letu GV be the unique solution 
of (a)-(h). Then it holds that u < w a.e. in LI for all w G Kq H Ki satisfying Aw G LP(LI). 


Proof. We set N := {x G LI: w(x) < u(a:)} and v := u A w. Since w G Kq, by Lemma 13.41 and 
da, V satisfies 


V G Kq, 



a.e. in Ll\N, 
a.e. in N, 


Vw 


Vu a.e. in Ll\N, 
Vw a.e. in N. 


Substituting v into the variational inequality (b) of Proposition 13.11 we have 



f){w — u) dx. 


0 < {Au — f,v — u)v 
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Since Au — / > 0 and w — u < 0 a.e. in N, it follows that Au = f a.e. in N. Here we note that 


V — u GV, V — u = 


0 


a.e. in n\N, 


V(n — m) = 


0 a.e. in n\N, 

\7{w — u) a.e. in N. 


w — u < 0 a.e. in N, 

From the fact that u — v = {u — w)+ and Aw G we obtain 

0 < a{u — v,u — v) = a{u — w,u — v) = {A{u — w),u — v)v 

= / [Au — Aw){u — w) dx = / [f — Aw)[u — w) Ax. 

Jn Jn 

Since Aw > f (by w G Ki) and u > w a.e. in N, we conclude that a(u — v,u — v) =0, whence 
u = V (hence u < w) a.e. in H. □ 


The lemma above will also play a crucial role for identifying the limit of each solution u = u[x, t) 
for (01)-O as t —oo in a proof of Theorem 12.91 
We are now ready to prove Theorem 13.91 

Proof of Theorem \S.!A By assumption, we find that C and Kf C K\. Moreover, U 2 
belongs to both and Kl := {v G V: Av > fi in H'}, and Au 2 G L^(H) as well. By Lemma 
13.101 we conclude that ui < U 2 a.e. in B. □ 


4. Reduction to an evolution equation and the uniqueness of solution 


In this section, we first reduce the problem (111)-(I7|) to the Cauchy problem for a nonlinear 
evolution equation in L^[Gt) with the aid of convex analysis. Then we shall prove Theorem 12.41 
on the uniqueness of solution. 

Let us begin with reformulating m as a parabolic inclusion with a multivalued nonlinear 
operator acting on the time derivative of u[x, t). Let a : M ^ 2® be given by 


a{s) 


{ 0 } 

(- 00 , 0 ] 


if s > 0, 
if s = 0 


(31) 


with the domain D[a) = [0,oo). Then s + a[s) is the (multi-valued) inverse mapping of the 
function (s)+, and it can be also represented by 


a{s) = 9/[ . >o](s) for s > 0, 

where /[ .>o] denotes the indicator function over the set [s > 0] := {s G K.: s > 0} and d means 
the subdifferential in the sense of convex analysis (see, e.g., [TS] and also (1551) below with iL = R). 
Then (0]) can be reformulated as a doubly nonlinear-type PDF, 


dtu a[dtu) 9 Am + f in Q. 


(32) 


We next reduce the PDF (1321) to an evolution equation. To this end, define a functional 
(j) : LF'[Vl) [0, oo] by 

|-|-oo ifMGL^(B)\B 

with the effective domain D[(j)) := {v G : (j)[v) < -|-oo} = V. Then we observe that: 

Lemma 4.1. The functional (f> is convex and lower semicontinuous in LF'[^1). In particular, if 
is satisfied, then the subdifferential operator d(j> of ft [in L^(B)) is characterized as 

D[d4>) = X nV, d4>[v) = —Av for v G X flV, 

where A stands for the Laplace operator from D[A) = X CiV into L^[Lt) as in Provosition AS. fii 
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Here let us recall the definition of the subdifferential operator dip : H ^ H oi a proper, lower 
semicontinuous and convex functional ip defined on a Hilbert space H, 

dip{u) := {^ G H-. p{v) — ip{u) > {f,v — u)h for all v G D{ip)} for u G D{ip), (33) 

where (•, ■)h stands for the inner product in H and D{ip) := {w G H: p{w) < +oo}, with domain 
D{dp) := {w G D{p): dp{w) 7 ^ 0}. It is well known that dp is a (possibly multivalued) maximal 
monotone operator in H (see, e.g., [H] for more details). 

Proof of Lemma \4.1\ We note that the restriction cfo := (t>\v of (j> onto V is Frechet differentiable 
and the derivative (j)^ of (fo satisfies 

{(f)Q{u), z)v = [ Vit-Vzda; for all z G V. (34) 

JQ 

Now, let u G D{d(j)) and ^ G d(j){u) C L^(H). From the definition of subdifferentials, we find 
that d(j){u) C d(j)o{u) = {(/<o(u)} for all u G D{d4>) C D{(p) = V. Hence = <('o(u), i.e., 
d(j){u) = {4>q{u)} and 4>'q{u) G i^(fl); here and henceforth, we simply write d(f){u) = 4 >q{u). It 
follows that 

Aiu = u + = u + G Lf{Pl). 

Moreover, by 03 along with Proposition 13.61 we deduce that u — {u - 1 - ^) G X Pi V and that 
u + f = Aiu = —Am + M. Therefore we deduce that d(j){u) = (j>Q{u) = —Am and D{d(f)) C X (IV. 
On the other hand, it is clear that X DV C Dl^df)), and hence, D^dcf) = X nV. □ 


Therefore the initial-boundary value problem for (I32|) equipped with (O-© can be rewritten 
as the Cauchy problem for an evolution equation in Lf{Vl) of u{t) := u{-,t)i 

dtu{t) + dl[ .>o]{dtu(t)) + d(f{u(t)) 9 f(t) in T^(0), 0 <t <T, m(0) = Mq, (35) 

where f{t) := f{-,t) and df .> 0 ] denotes the subdifferential operator in L^( 0 ) of the functional 
/[ .> 0 ] : Lf{VL) [0,oo] defined by 




0 if M > 0 a.e. in O, 
00 otherwise 


for V G 


We note that dl[ .>o](^^) = <a(^^(')) ^ where a(-) is a multivalued function given by 

(1^ . and D{dl[ . > 0 ]) = {m S : v >0 a.e. in 0} (see, e.g., [H]). 

Here and henceforth, for simplicity, we use the same notation /[ . > 0 ] for the indicator function 
over [0, -foo) defined on R as well as for that over the set {v G L^(H): m > 0 a.e. in fl} defined 
on L^(H), unless any confusion may arise. Moreover, the subdifferential operators of the both 
indicator functions are also denoted by dl[ . > 0 ]. 

Strong solutions of (1551) are defined as follows: 


Definition 4.2 (Strong solution of (1551) '). For given f G L‘^{Qt) o-nd uq G L^(H), a function 
u G (^([O, T]; T^(n)) is called a strong solution of (I11)-(I7]) on [0,T], if the following conditions 
are satisfied: 

• uGW'^'^{0,T;L'^{n))nL^{0,T;Xr\V); 

• It holds that 

dtuff) + dl[ . >o](i9tM(t)) + d(j){u[f)) 9 f{t) in for a.e. t G (0,T); (36) 

• u(0) = Uq, 

where the functionals /[ .> 0 ] and 4> on L^(fl) are defined as above. 

Proposition 4.3 (Equivalence of solutions). The notion of strong solutions for (I35II is equivalent 
to that for (fT |) -([7 1) defined by Definition \2.2[ 


14 


GORO AKAGI AND MASATO KIMURA 


Proof. Since a{s) is the inverse mapping of s i-A (s)+, one observes that (HD is equivalent to (1321) 
at each (x, t) G Q. Moreover, due to Lemma Hdl for each strong solution u of (1551) . u{x, t) satisfies 
(15^ a.e. in n X (0, T). Conversely, let it be a strong solution of ([T|) in the sense of Definition 12.21 
Then from the regularity condition (ii) of Definition 12.21 the right-hand-side of the inclusion 

a{dtu) B Alt + f — dtu 

belongs to L'^{H) for a.e. t G (0,T). Moreover, recalling that 

dl[ ■>o]{v) = {^(•) G L^(n): ^(x) G a{v{x)) for a.e. a; G fl} 

(see above), the evolution equation (l36l) holds in for a.e. t G (0,T). □ 

We next provide a chain-rule for the function t i—>■ 4>{u(t)), which is derived from a standard 
theory on subdifferential calculus and which will be used frequently to derive energy estimates in 
later sections. 

Lemma 4.4. We suppose that u G T; L^(D)) D L^(0,T; X (d F). Then we have'. 

(i) the function 

t i-B (f{u{t)) = I- [ \Vu{t,x)\'^ dx 

2 Jq 

belongs to W^’^(0,T); 

(ii) for a.e. t G (0,T), it holds that 

^ ^ = - J^dtuAudx, 

where denotes the inner product of 

(hi) itGC'([0,r];F). 

Proof. Thanks to Lemma 3.3 of [15] . the assertions (i) and (ii) follow immediately. Concerning 

(iii), since it belongs to L°“(0,T;F) and (^([O, T]; L^(fl)), by exploiting Lemma 8.1 of [32], one 
finds that it is continuous on [0,T] with respect to the weak topology of V. On the other hand, 
t I—> ||it(t)||v is continuous on [0,T] by (i). Therefore from the uniform convexity of || ■ ||v, we 
deduce that t >->• u{t) is continuous on [0, T] with respect to the strong topology of V. □ 

Before proceeding to a proof of Theorem 12.41 let us note that 

|a+— 6+p < |a+— 6+||a — 6| = (a+— 6+)(a — 6) for all a, 6 G M, (37) 

since the function s i—= s VO is nondecreasing and non-expansive, that is, |a+ — &+| < |a — 6|. 
Now, we are ready to prove Theorem 12.41 


Proof of Theorem \2.4\ For each i = 1, 2, let Ui be a strong solution of (|4|)-([7|) with uq = mq.z S F 
and f = fi G L^(0,r;L^(D)) and set u = ui — it 2 . Then due to Lemma l4^ we have 


l_d 

2^ 


iVup dx = 


l37t 

< 


which implies that 


I dtu Alt dx 


n 

f [(Am -k /i)+ - (Am + f2)+] [(Am -b /i) - (Am + /2) - /i + /2] dx 


- f |(Aui + /i)+ - 

j r2 

- l/l - f2\^dx 


(Am 2 + /2)+|^ dx 
for a.e. t G (0, T), 


/ |5tm - 9tM2p dx-b ^ / jVm - VmP dx < [ |/i-/ 2 pdx 
Jn “ Jn Jn 


for a.e. t G (0, T). 
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Integrate both sides with respect to t to obtain 

[ \\dtui{t) - dtU 2 {t)\\l 2 ,mdt+ sup \\Vui{t) 
Jo te[o,T] 




< 2 ^||Vuo,i - Vuo,211^2(0) + J 

In particular, if Mpp = uo ,2 and fi = / 2 , then ui coincides with 
solution of ©-(El) is unique. 


- f2(t)\\l2(^a)d?j . (38) 

U 2 a.e. in Qt- Consequently, the 

□ 


Corollary 4.5 (Continuous dependence of solutions on data). For each T > 0 and i = 1,2, let 
Ui be the strong solution of 0-0 on [0,T] with uq = u^ i G V and f = fi G L^{Qt)- Then 
(1551) holds true. 


5. Existence of solutions and comparison principle 


In this section, we shall prove Theorem l2.6l on the existence of solutions for (01)-©. Let T > 0 
be fixed. We denote by r a division {to,ti, ... ,tm} of the interval [0,T] given by 

0 = to < h < ... < tm = T, Tk :=tk - tk-i for fc = 1 , ..., m, | t | := max t ^. 

We shall construct Uk G X DV (for k = 1,2,..., m), which is an approximation of u{tk) for a 
solution u of m by the backward-Euler scheme 


Uk Uk — l , . I f ^ • r> 

- = (Auk + fk)+ a.e. m i2, 


(39) 


where fk G Lf (O) is given by 

fk-.= - r f{;s)ds. 

For given uq G E, we shall inductively define Uk G V for /c = 1, 2,..., m as a (global) minimizer 
of the functional 

Jk{v) f \v\'^ dx [ \Xv\'^ dx - + fk,v\ for u G E (40) 

Jo, ^ Jfl \ Tk / y 


subject to 


V G Kq := {v gV ■. V > Uk-i a.e. in f2}. 


(41) 


Remark 5.1 (Derivation of the discretized problems). The minimization problems with con¬ 
straints stated above can be also derived from a discretization of the evolution equation (1351) . 
which is equivalent to (|4]) (see Remark |2.31 and Section |4]). A natural time-discretization of (l35l) 
may be given as 

-—- ^»ud[.>o]l -- j-Auk5fk in E (42) 

(here dv stands for the subdifferential of the functional /[ . >o] restricted onto E), which is an 
Euler-Lagrange equation for the functional 

-E^fc(T) := ^ / |upda;-k/[ .>0] ^ ^ ^ 

Indeed, since Ek is coercive, lower semicontinuous and convex in E, Ek admits a global minimizer 
Uk over E, and moreover, Uk solves (j42ll in V'. Here we note that the minimization of Ek over E 
is equivalent to that of Jk over Kq from the fact that 


> 0 ] 




0 if ?; > Uk-i a.e. in H 
oo otherwise 


for V G L/^iTl). 
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Applying the regularity theory established in Section |3l one can actually obtain the unique 
minimizer Uk of Jk over for each k. More precisely, we obtain the following theorem, where 
we set 


gk '■= 


Tk 


^Uk - fk- 


Lemma 5.2 (Existence and regularity of minimizers). For given uq G V and each fc = 1, 2,..., m, 
there exists a unique element Uk G Kq which minimizes dMl) subject to SB)- Moreover, for each 
k = 1,2,... ,m, the minimizer Uk belongs to X and fulfills dMl), that is, 

Uk — Uk-i >0 a.e. in fl, (43) 

gk > 0 a.e. in O, (44) 

{gk,Uk-Uk-i)v = Q, (45) 

Furthermore, one has 

{gk,v - Uk)v > 0 for all v £ K^, (46) 

(Sfc +/fe + Aufc-i,u - Mfe)y > 0 for all v G Ki, (47) 


where A is the Laplace operator from X nV into Lf{Ll) {see Section^ and the set is given 
by 

K^-.= [vGV-. ^ ~ -Av-fk>0 in V' 

[_ '^k 

In addition, suppose that the conditions (jl8|) and (|2^ hold. Ifuo £ and {/fcj^i C L^{Lt), 

then {ukj^^i Cl XP and it holds that 

0 < < {—Auk-i — fk)^ OL.e. in ft for each k = 1,... ,m. (48) 


Proof. Let us start with fc = 1. By setting a = l/r^, > 0 (i.e., Au = A^u = u/xk — Au), 
f = fk + Uk-i/xk £ L‘^{Cl) and ^ = Uk-i C X PV, one can write Kq = Kq and Jk{v) = J{v) 
for u £ E with Kq and J{v) defined by (IT4l) and (ITSl) along with (fT^ . Hence, one can apply 
ProDOsition IS . 1 1 and Theorem l3.2l to the minimization problem of Jk over Kq. Then the minimizer 
Uk £ Kq of Jk over Kq uniquely exists, and furthermore, (l4^ -(IT7 )l follow immediately from the 
fact Uk £ Ki, (b), (c) and (d) of Proposition 13.11 respectively. Moreover, by virtue of (fTOll and 
Corollary 13.81 one can deduce that Uk £ X. Repeating the argument above for fc = 2, 3,..., m, 
we can inductively obtain Uk £ Kq n X satisfying (P l) - (|Tf)) for each fc = 2, 3,..., m. 

Finally, if uq £ X^ nV, fk C LP{Cl) and (E51) is satisfied for p satisfying (IT^ . by Corollary 13.81 
we can assure that Uk £ X^ n Kq and 

fk + < — - Auk < (fk + V f'Hflzl _ Auk-i) for a.e. x £ H, 

Tk Tk \ Tk J \ Tk ) 

which is equivalent to (l48l) . □ 


Proof of Theorem \2.6[ Let us define the piecewise linear interpolant Ur £ W^’°°{Q,T',X n V) of 
{uk} and the piecewise constant interpolants Ur £ L°°{0,T-,X n V) and fr £ L°“(0,T; L^(0)) of 
{uk} and {fk}, respectively, by 

Ur{t) := Uk-i + - —^^(ufc-Ufc-i) for tG[tk-i,tk] and fc=l,...,m, 

'^k 

Ur{t):=Uk, frit) ■■= fk for tG{tk-i,tk] and fc=l,...,m. 
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By summing up (I45p for /c = 1,..., £ with an arbitrary natural number £ < m, we have 
2 „ 


E 

fe=i 


Tfe 


Uk — Uk-1 


Tk 


+ - / \Vwfdx-- [ iVwol^dx 


L2(n) 


In 


In 


- E'^'= ( 

k=l 

which implies 


'^k Uk—1 
Tk 


L 2 (a) 


£ £ 


Tk 


k=l 


k=l 


Uk — Uk-l 


Tk 


(49) 


L 2 (n) 


J l| 5 tMT(s)||^ 2 (f 2 ) ds + IIV'Ui-(t)||i2(f2) < II V'Uo||i 2 (Q) + J ||/T(s)|y 2 (Q) ds for all tG[0,T]. 
Hence, we obtain 


ll^tU,- 11 ^ 2(0 T’. ^2 (Qp + sup II VuT-(^)llL 2 (n) d” ll^'*^'r( 0 llL 2 (f 2 ) 

' ^ " te[0,T] ^ ' te[0.T] ' ' 


< c 


2 + II f lp 

^ \\J'^\\L^{0,T-L^{n)) 


L 2 (n) 


)■ 


(50) 

(51) 


Now, let us take a limit as m ^ oo such that |r| —>■ 0 and note that 

Jr ^ f strongly in L^( 0 ,r;L 2 (H)). 

In particular, {/,-} is bounded in i^(0, T; L^(H)). Indeed, one can verify that 

ll/T||L2(0.T;L2(n)) < ||/||L2(o_T;L2(n))- 

From the uniform estimate (ISO)) , one can take a function u G IF^’^(0, T ; i^(H)) n L°°{0, T ; V) (in 
particular, u G C{[0,T]; L^{V,)) as well) such that, up to a (non-relabeled) subsequence, 

Ur ^ u weakly in kF^’^(0,T;L^(f2)), (52) 

weakly star in L°°(0,T;H), (53) 

strongly in C([0, T]; ^^(H)), (54) 

Ur ^ u weakly star in L°°(0,T;H), (55) 

Ur{T) —>■ u{T) weakly in V. (56) 

Here, the weak and weak star convergence of Ur and Ur immediately follow from the uniform 
estimate (1501) . Moreover, we also note that Ur and Ur possess a common limit function. Indeed, 
by a simple calculation, we observe that 

tk t 


lillr(^) 7/7-(t) II L2(n) — 


< 


B9t 


Tk 

Uk — Uk-1 


\\uk Ufc-l|li 2 (f 2 ) 


Tk 


Tk 


L 2 ( 0 ) 


which yields that 


< for all tG ( 4 - 1 , 4 ], fc = l, 2 ,..., 


sup ||M 7 -(t) - Mr( 0 llL 4 n) < C|r|^/^ -)> 0 . 
te[o,T] 


Thus Ur and Ur (weakly) converge to a common limit function. Furthermore, since V is compactly 
embedded in L^(H), due to Ascoli’s compactness lemma along with (1501) . we obtain the strong 
convergence (iMl) . Since Ur{T) = Um is bounded in V by (l4^ . one can also derive (l56l) from (l54)) . 
We further observe that m(0) = uq. 

We next estimate Aur in 4^(0, T; T^(H)) by using (l4^ and the assumption (fTTl) . We first 
rewrite (H51) as 


Uk - Uk-1 
Tk 


+ fk < -Auk < - 


'^k Uk—1 
Tk 


+ /fe W — 


Uk — Uk-1 
Tk 


- Auk-1 


a.e. in H. (57) 
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Since {uk — Uk-i)/Tk > 0 a.e. in il hy Uk € Kq, we observe by (fTTI) that 

(The right-hand side of (1571) ') < fk^ (—< /* V {—Auk-i) a.e. in fl, 
which also iteratively implies that 


-Auk < f* V {-Auk-i) 

< /* V if* V i-Auk-2)) 

= f* i—Auk- 2 ) <■•■</* V (—Alto) a.e. in fl. 


Thus we obtain 


Uk Uk—1 
Tk 


+ fk< -Auk < f* y (-Alto) a.e. in Q, 


which yields that 

Il^'afc|li2(n) - 2 ( Il/*lli2(n) + ll^''^o|li2(n) 


Uk - Itfe-l 


Tk 


1,2(0) 


ll/fellL2(0) 


for k = 1,2,... ,m. Hence we deduce that 

J II^Ur(0lli2(o) < 2T ^||/*||i2(o) + ll^wo|||2(o)) 

+ 2y^ \\dtUrit)\\l2(^a) dt + 2 \\frit)\\l^^^^ dt < C (58) 

by using (ISU)) and (ISTl) . 

Exploiting Proposition l3.6l with (ITOl) . we see that (/ — A) G Isom(AnF, T^(fl)), which together 
with (|5^ gives 

^ \\Ur{t)\\x dt < C dt < C. 

Therefore we have, up to a (non-relabeled)subsequence, 

Ur ^ u weakly in L^(0, T; A), 

Aut-—A lt weakly in L^(0, T; T^(H)), 

which particularly implies it(t) G D{A) = X nV for a.e. t G (0,T). Therefore the piecewise 
constant interpolant gr of {gk} defined by 

grit) := gk ^ - Aitfc - fk for t G (4-i,4] 

'^k 

converges to 

dtu — Alt — f =■ g (59) 

weakly in T^(0, T ; T^(H)). 

It remains to prove that it solves (|4]) for a.e. ix,t) G Qt- To this end, we recall the evolution 
equation (1551) equivalent to dH). Then it suffices to check that 

dtu>0 a.e. in Qt and — g(t) G i9/[ . > 0 ] (9tit(t)) for a.e. t G (0, T). 

To this end, we employ the so-called Minty’s trick for maximal monotone operators, since dl[ ■ > 0 ] 
is maximal monotone in L^(H). 


Proposition 5.3 (Demiclosedness of maximal monotone operators (see, e.g., [T51 [T51 [TU])). Let 
A : H ^ H be a ipossibly multivalued) maximal monotone operator defined on a Hilbert space H 
equipped with a inner product (•, •)//. Let [un,^n\ be in the graph of A such that Un ^ u weakly 
in H and £,n ^ f weakly in H. Suppose that 

limsup(^„,M„)// < i^,u)H- 

n—¥-\-oo 
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Then [m, belongs to the graph of A, and moreover, it holds that 

lim {^n,U„)H = 

n—^+oo 


Note that {uk — Uk-i)/Tk > 0 a.e. in For an arbitrary w € D{I^ • >o]) = {'P e L^{n): V > 
0 a.e. in substitute v = WTk + Uk-i £ Kq to (HSl) . Then we see that 


Bell 

0 > {-gk, V - Uk)v = Tk ( -fffe, w - 


Uk — Uk-1 
Tk 


L^(Q) 


which together with the arbitrariness oi w £ D{I[ . >□]) and the definition of /[ . >o] implies that 

-gk G dll . > 0 ] ~ ^ ■ 


Tk 

Moreover, for k = 1,2,... ,m, we find by (ITil) that 

2 




Uk — Uk-l 


Tk 


< - 


L^(Q) 


'^k '^k — 1 


Tk 


(j){uk) - cfiiuk-l) 


z,2(n) 


Tk 


fk 


Uk — Uk-l 
Tk 


L2(n) 


which leads us to get 

pT 


/o 


(-5rW,atUr(t))i2(o) dt<-f ||l9tUr(t)||i2(f2) dt - (T)) + </)(uo) 

Jo 


+ {fT{t),dtUrit))^^,^. dt. 


L2(n) 

^ \j 

Taking a limsup as | t | —^ 0 in both sides, exploiting the weak lower semicontinuity of norms and 
the functional and recalling Lemma [4.41 we conclude that 

limsup / [-gr{t),dtUr{t))^ 2 i^) dt<- ||9tM(t)|li2(n) dt - (j){u(T)) + (j){uo) 

ItMo Jo Jo 


<f^ 


/ {f{t),dtu{t))^2^a) di 
0 

{-dtu{t) + Au{t) + f{t), 5tu(t))^2(o) dt 

■T 

i-g{t),dtu{t))^2(^a) dL 


(60) 


Consequently, by virtue of the (weak) closedness of maximal monotone operators (see Propo¬ 
sition [O] above), it follows that dtu{t) G D{dli .>o]), i-e., dtu{t) > 0 a.e. in ft, and —g{t) G 
dl[ ■>o](^t^(0) foi' S'-®- ^ € (.d,T). Therefore u solves (1551) . and hence, u is a strong solution of 
O-Gl)- Thus Theorem [5751 has been proved. □ 


Remark 5.4. To prove that m is a strong solution of 0-0, it is possible to show the conditions 
(V1)-(V6) of Theorem l7.1l in instead of the last argument of the proof of Theorem l2.61 Actually, 
(V2) and (V3) directly follow from (1431) and (|551) by taking limit of |r| —>■ 0, respectively. The 
condition (V4) follows from the estimate (160L since the left-hand side of ((5D1) is zero and the 
right-hand side is non-positive. 


Due to Theorem [211 the limit of {ur} and {ur} is unique, whence they converge along the 
full sequence. 

Corollary 5.5. Seguenees {u,-} and {mt} converge to the unigue solution u of 0-0 as |r| —;■ 
0 +. 


We next prove Theorem 12.81 
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Proof of Theorem \2.8[ Let and be strong solutions of (|4]) with uq = Uq and / = /* for 
i = 1,2, respectively. By the uniqueness of solutions (see Theorem 12.41) and the construction 
of solutions discussed so far, one can take discretized solutions {u].} for i = 1,2 such that the 
piecewise linear interpolant u\. of { m ^} converges to u® strongly in C([0, T]; L^(0)) as | t | —0, 
and they solve the variational inequalities 

u\. S Kq{u\_i) •.= {v & V V > Mfe_i a.e. in fl}, 

ak{ul,v-ul)> f {fl + ul_j^/Tk) {v-uDdx lor all v G Ko{ul_i), 

Jn 

where afe(-, •) stands for the bilinear form given by 

ak{u,v) = I Vu-VudxH - f uvdx for u,vGV. 

Jq Tk Jq 

By iteratively applying the comparison theorem for elliptic variational inequalities (see Theorem 
[HI), from the fact that a.e. in Qt and < Uq a.e. in 17, one can deduce that 

Mfe < a.e. in Qt for all fc = 1, 2,..., m, 

which also implies u].{t) < u^(t) a.e. in 17 for all t G (0,T). Then passing to the limit as |r| —0, 
we conclude that < v? a.e. in Qt- D 


6. Long-time behavior of solutions 


This section is devoted to proving Theorem [231 Let us begin with deriving a uniform estimate 
for u{t) for i > 0. To do so, recall the construction of the unique solution u = u{x,t) of (|ll)-(l2l) 
performed in the proof of Theorem l2.6l and particularly note that 


Uk > Uk-i and gk ■= 

It follows that 

hk '■= Qk d- fk = 

We also recall the estimate 


Uk — Uk-l 


Tk 

Uk — Uk-l 


— Auk — fk A 0 a-e. in 17. 


— Auk > —Auk a.e. in 17. 


Tk 

which gives 
fk <hk < {-Auk-i) V fk a.e. in 17. 

Therefore by (H3) we find that 

163 } 111 } 

fk < hk < {-Auk-i)V fk 

lea 

< hk-iVf* 

< {hk-2 V n V r 

lea 

= hk -2 V /*<•••< hi V /* < (-Alto) V f* a.e. in 17, 
which together with the assumption that uo G X llV gives 

\\hk\\L‘^{Q) < ||AMo|U2(n) -f ||/*||l2(q) -I- ||/fc||L2(Q) 

< ||AMo||L2(n) + ||/*||L2(n) -I- ||/||Loo(o_oo;L2(n)). 

Here we used the fact that ||/fc||L 2 (Q) < ||/||L~(o,oo;T 2 (n)) for all k. Moreover, set 
hr{t) := dtUrit) - Aur{t) = hk for t G {tk-i,tk]. 

Recalling the convergence of approximate solutions obtained in the proof of Theorem 
observe that 

hr —>■ dtu — Alt =: h weakly in L^(0, T ; L^(17)), 


(61) 

(62) 

(63) 


(64) 


we 




UNIDIRECTIONAL EVOLUTION EQUATIONS OF DIFFUSION TYPE 


21 


On the other hand, since {hr} is bounded in L°°(0, T; L^(f2)) by (IMl) . we assure, up to a (non- 
relabeled) subsequence, that 

hr ^ h weakly star in L°°(0, T; L^(f2)) 

as |t| 0. Moreover, from the lower semicontinuity of the L°°-norm in the weak star topology, 

we have, by (IMl) . 

ll^llL“(0,T;L2(n)) < li^nf < ||AMo||L2(n) + ||/*||L 2 (n) + ||/||loo(o.oo;L 2 ( 0 )) 

for each T > 0. Since the bound is independent of T > 0, one can derive that 

ll^llL“(0.oo;L2(n)) < ||Auo||l 2 (q) + ||/*||l 2 (q) + ||/||Loo(0.oo;L2(n))- (65) 

Note that (C)'L')z, 2 (o) = 0 for all v G L^(0) with w > 0 a.e. in O and ^ G dl[ .>o](^^)- Thus 
testing (I55|) by dtu{t), we have 

119*^(01112(0) + ~\\ ynmUn ) = (/(0,9*u(0)l2(o) 

= (/(O - /oo,9tu(0)L2(o) + ^(/oo,w(0)L2(n) 

< ^l|9tw(0lli2(O) + ^11/(0 -/oo||i2(o) + ^(/oo,u(0)L2(f2). 
Define an energy functional E on V hy 

^^(^) := ^l|Vw||i 2 (n) - (/oo,?^)L 2 (n) for v gV. 


Then one has 

^I|9tw(0lli2(o) + ^^^(w(0) < ^11/(0 -/oo||i 2 (o) for a.e. t > 0, 
which implies the non-increase of the function 

t ^ E{u{t)) - ^ II/(t)-/oo 11^2(0) dr for t > 0. 
Moreover, by using the Poincare inequality (due to (HI)), we have 

Eiv) > Jl|Vu||i.(^) - C||/oo||i 2 (a) for all u G H. 

Thus integrating (l66)) over (0, s) and using (H2) and (l67l) we obtain 

pCO 

/ l|9tM(0lli2(n)dt < C, 

^0 

SUp||V-li(t)||i2(Q) < C, 
t>o 

which also yields 


sup IIAM(0||y' < C. 

t>Q 


( 66 ) 


(67) 


( 68 ) 

(69) 


(70) 


Moreover, by virtue of (USD, we have 

I|5||l“(0,oo;L2(O)) < ||^||l“(0,oo;L2(O)) + || / || (0,oo;L2 (Q)) < M (71) 

for some constant M. Here we used that g{t) = h{t) — f{t). 

Let I C (0,oo) be the set of all t > 0 for which (1551) holds true and ||g(0llL2(o) is bounded by 
M as in (ED). Then the set (0, oo) \ / has zero Lebesgue measure. Recalling by (1551) and (H2) 
that 

^ (l|9tu(i)||i2(o) + Wfit) - /oo||i 2 (o)) dt < oo, 
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one can take a sequence £ [n, n + 1] n / such that 

dtu{sn) 0 strongly in (72) 

f{sn) /oo strongly in (73) 

as n —>■ oo. 

Moreover, by using the preceding uniform (in t) estimates and using the compact embedding 
V ^ 77^(0), we deduce, up to a (non-relabeled) subsequence, that 


u{sn) —t -z weakly in F, 

(74) 

strongly in 

(75) 

Au{sn) —t —Az weakly in V', 

(76) 

—g{t) —>■ / weakly in Lf{Vl) 

(77) 


with some z € V and ^ From the demiclosedness of dl[ . >o] in L^(fi) and the fact that 

—g{t) £ 9/[ .>o](9tM(t)) for a.e. t £ (0,oo), it follows that ^ £ dl^ ■>o](0); that is, ^ < 0 a.e. in 
n by dl[ .>o](0) = (—00,0]. Moreover, by (I55|) and (1751) . we get ^ — Az = foe, which leads us to 
/oo + Az = / < 0 a.e. in il. Furthermore, by (ITU)) along with Proposition 13.61 the limit z belongs 
to X, since z — Az = 2 — / + /oo £ i^(f2). 

Therefore we derive that 

l|Vu(sn)||i2(n) = (-Au(s„),u(s„))^2(n) 

= { — dtu{Sn),u{Sn)) : u{Sn)) ]^2 + (/(Sn) , M(Sra))^2(Q) 

^ (“? + /oo, ^)j^2(f2) = II Vz||^2(f2). 

From the uniform convexity of F, it holds that 

u{sn) —t strongly in V. (78) 

We shall next verify the convergence of the solution u{t) to the same limit z as t —>• oo, that is, 
u(tn) —>• z for any sequence —)• oo and the limit z is independent of the choice of the sequence 
{tn). Subtracting the stationary equation 

dl[ ■ >0] (0) - Az 9 /oo 

from the evolution equation (1351) . we see that 

dtuft) + df . > 0 ] idtu{t)) - 5/[ . > 0 ] (0) - A {u{t) - z) 9 f{t) - foe- 
Test it by dtu{t) to get 

\\\dtu{t)\\l2(^) + ^^l|V(u(t) - Z)||i2(0) < \\\f{t) - /oo||i2(0). 

Integrate both sides over {sn,T) for t > Sn- Then it follows from (1781) that 
1 11 

- sup ||V(M(r) - z)||i 2 (o) < -||V(m(s„) - z)||i 2 (n) + - / ||/(t) - /oo|li 2 (n) dt 

Z T>s„ z L J 


Thus u{t) converges to the limit z strongly in F as t ^ oo. This completes the proof of the first 
half of the assertion. 

We next prove the second half of the assertion. In addition, assume that f{x,t) < foe{x) for 
a.e. X € Q and let z £ A D F be the unique solution of the variational inequality (Vl)(uo, /oo). 
Then by Proposition 13.II and Theorem 13.21 for A = Aa- with ct = 0, z satisfies — Az > foe a.e. in 
n, and moreover, we deduce that U{x,t) := z{x) becomes a strong solution of (H)) by observing 
that 

dtU = 0 and A17(a:, t) +/(a:, t) < Az(a;) +/oo(a:) < 0 a.e. in Q. 
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Hence by the comparison principle for the evolutionary problem (|4]) (see Theorem l2.8|) . we assure 
that u{x,t) < z{x) for a.e. {x,t) G Q. Letting t > oo and recalling (TTSl) . we obtain 

z{x) < z{x) for a.e. a: G H. 

On the other hand, since z belongs to X and satisfies z > uq and —Az > foo in V, 
applying the comparison theorem for variational inequalities of obstacle type (see Lemma 13.101) 
to (VI)(uo; /oo), we assure that z < z a.e. in O. Consequently, we conclude that z = z a.e. in O. 
Thus we have proved the second half of the assertion of Theorem 12.91 □ 


7. Other equivalent formulations 

In this section, we discuss other formulations of solutions for 0-0 equivalent to those defined 
by Definition 12.21 Let us start with a complementarity form of strong solutions. 

Theorem 7.1. Let f G L‘^{Qt) and uo G L^(0). Then u is a strong solution of the problem 
0-0 on [0,T], if and only if the following six conditions are satisfied: 

(VI) u G w^'^{o,T- L^{n)) n L^{o,T- H^{n) n H^n)), 

(V2) dtu >0 a.e. in Qt, 

(V3) dtu — Au — / > 0 a.e. in Qt, 

(V4) {dtU — Au — f) dtU = 0 a.e. in Qt, 

(V5) = 0 'H^~^-a.e. on Ln and u = 0 'HT~^-a.e. on Td for a.e. t G {0,T), 

(V6) m( 0, •) = Mo- 


Proof. If u satisfies (VI) (or (i) of Definition [50, one can define the following measurable subsets 
of Qt- 


Qo ■= {{x, t) gQt: dtU ^ {Au + /)+}, 

Qi ■= {(a:, t) &Qt: dtU = {Au + /)+ > 0}, 

Q 2 ■■= {{x, t) gQt: dtU = (Am + /)+ = 0}, 

which are disjoint and satisfy Qt = Qo U Qi U (52- 

Let u satisfy (i)-(iii) of Definition 12.21 Conditions (VI), (V5) and (V6) follow immediately. 
From (ii) of Definition 12.21 it follows that 


7^"+i(Qo) = 0, 

and moreover, by definition, 

dtu > 0, dtu — Au — / = 0 a.e. in Qi, 
dtu = 0, dtu — Au — / > 0 a.e. in Q 2 . 


Hence (V2), (V3) and (V4) follows. Consequently, every strong solution u of 0-0 in the sense 
of Definition 12.21 satisfies all the conditions (V1)-(V6). 

Conversely, let u satisfy (V1)-(V6). Conditions (i) and (iii) of Definition [20 follow from (VI), 
(V5) and (V6). Let us next show that 'H'^~^^{Qo) = 0, which is equivalent to the condition (ii) of 
Definition 12.21 Define 


Q* := {{x,t) G Qt'. dtU > 0 and dtU — Au — f > 0 at {x,t)}. 

Then it holds that 'H^~^^{Qt \ Q*) = 0 by (V2) and (V3). 

We claim that 


dtu > 0, dtu — Au — / > 0 at each {x, t) G Qo n Q*. 


(79) 
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Indeed, by the definitions of Qo and Q*, u satisfies the following conditions at each (x, t) G QoHQ*: 


dtu ^ (AM-h/)+, 

(80) 

dtu > 0, 

(81) 

dtu — Au — / > 0. 

(82) 


If dtu = 0 at some (xo,to) G Qo^Q*, then Au + f > 0 at (a;o,to) by (l80l) and it contradicts (l82|) . 
Hence, we obtain dtU > 0 at each point of Qo n Q* by (ISTI) . Similarly, if dtU — Au — / = 0 at 
some (a;o,?/o) € Qo<^Q*, then 0 < dtU = Au + f = (Au + /)+ at (xo,yo), which contradicts (|5(I)) . 
Thus, we obtain dtU — Au — / > 0 in Qg H Q* by (1^ . 

Since Qt = Qo U Qi U (52 is a disjoint union and \Q*) = 0, we have 


0(V4) 


Qt 


{dtu — Au — /) dtu dx dt = / / {dtu — Au — /) dtu dx dt 

J Jqo 


{dtu — Au — f) dtu dx dt. 

iQonQ* 

By (17^ . we obtain H"+^((5o n(5*) = 0; otherwise the last integral is positive. Hence we conclude 
that 


H”+^(Qo) = H"+^(Qo n Q*) + \ Q*) = 0. 

This completes the proof. 


□ 


Finally, let us discuss a possible formulation of (HI) in the sense of viscosity solutions. Set 

F{x,t,Y) :=-{tvY + f{x,t))^ for x G H, t G (0,T), F G (83) 

where denotes the set of all symmetric n x n real matrices. Then o is also written as 

dtu{x, t) + F(x, t, D'^u{x, t)) = 0, 

where D‘^u{x, t) G R"y^ is the Hessian matrix of u. Since F is degenerate elliptic, one may apply 
the theory of viscosity solutions to prove the existence and uniqueness of viscosity solutions of 
(1831) under suitable assumptions for f{x, t) and the boundary condition. However, to the authors’ 
knowledge, no result on such a viscosity approach to has been obtained except for [li] . 
Moreover, the relation between the notion of viscosity solutions and that of strong solutions for 
(HI) is widely open. For further details of the theory of viscosity solutions, we refer the reader 
to m. EH, m and references therein. 
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